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Recently, two independent calculations have been presented of finite- mass ( "self-force" ) effects on 
the orbit of a point mass around a Schwarzschild black hole. While both computations are based on 
the standard mode-sum method, they differ in several technical aspects, which makes comparison 
between their results difficult — but also interesting. Barack and Sago [Phys. Rev. D 75, 064021 
. . , (2007)] invoke the notion of a self-accelerated motion in a background spacetime, and perform 

OO ■ a direct calculation of the local self-force in the Lorenz gauge (using numerical evolution of the 

' perturbation equations in the time domain); Detweiler [Phys. Rev. D 77, 124026 (2008)] describes 

' the motion in terms a geodesic orbit of a (smooth) perturbed spacetime, and calculates the metric 

perturbation in the Regge- Wheeler gauge (using frequency-domain numerical analysis). Here we 
^ ' establish a formal correspondence between the two analyses, and demonstrate the consistency of their 

O ' numerical results. Specifically, we compare the value of the conservative 0(/i) shift in (where jj, is 

, the particle's mass and is the Schwarzschild t component of the particle's four- velocity), suitably 

mapped between the two orbital descriptions and adjusted for gauge. We find that the two analyses 
^— »' ' yield the same value for this shift within mere fractional differences of 10~^-10~^ (depending on 

' , the orbital radius) — comparable with the estimated numerical error. 

U: 

O^, I. INTRODUCTION 

u ; 

Two recent works discussed the effects of gravitational self-force (GSF) on the motion of a point particle in a 
] circular orbit around a Schwarzschild black hole. Barack and Sago [1,] (hereafter BS) carried out a direct calculation 
Ps) . of both dissipative and conservative components of the local GSF in the Lorenz gauge. Detweiler [2] (hereafter SD) 
J> ' focused on quantifying a number of gauge-invariant orbital effects associated with the finite mass of the particle. 
Both analyses employ the mode-sum method [1, 0] , wherein the metric perturbation associated with the particle is 
first decomposed into multipole harmonics, and the contribution to the physical back-reaction force is then calculated 
, mode by mode. However, the two analyses use very difference languages in describing the motion of the particle. 
• ■ BS's description adheres to the original fundamental formulation of Mino, Sasaki and Tanaka 0] (also Quinn and 
' Wald (6| and the recent Gralla and Wald [?!), wherein the GSF is viewed as accelerating the particle with respect to 
a background geodesic. SD, instead, depicts the particle's orbit as a geodesic of a smooth, perturbed spacetime — an 
interpretation first put forward by Detweiler and Whiting Q. The two descriptions of the motion are known to be 
fundamentally equivalent [§] (see [sl] for a pedagogical review), but the different language they use makes it rather 
difficult to compare between the results of the respective calculations. 

The goal of this report is to establish a "common language" to facilitate comparison between the results of BS 
■ and SD, hence show that the two sets of results are in full agreement with each other. The motivation for such a 
comparison is three-fold: Firstly, it provides a reassuring confirmation for the equivalence between the two descriptions 
of the motion. Secondly, it demonstrates how the correspondence between the two descriptions is to be established in 
practice — a similar technique could then potentially be implemented in future calculations of the GSF. Thirdly, the 
comparison provides a good overall check of both calculations. It is indeed a good check, because the two analyses are 
highly independent: Not only do they invoke different interpretations of the perturbed motion, they also use different 
gauges for the metric perturbation, and apply very different numerical methods to evaluate it. 

The GSF itself is gauge dependent and, of course, so is the metric perturbation. To compare the results of two 
calculations held in different gauges, one naturally seeks to consider gauge-invariant quantities which are affected by 
the finite mass of the particle in a non-trivial way. A list of such quantities, for circular orbits in Schwarzschild, was 
introduced by Detweiler in [ll| (cf. SD). In discussing gauge-invariant GSF effects, it is useful to distinguish between 
effects which are dissipative and others which are conservative. Dissipative effects are manifest in a gradual, secular 
drift in the value of the intrinsic orbital parameters (such as the energy and angular momentum) due to the emission 
of gravitational waves. In the case of a circular equatorial orbit in Schwarzschild, such effects arise from the t and 
ifi components of the GSF (henceforth t,r,9,Lp denote the standard Schwarzschild coordinates), which are directly 
related to the rate of change of orbital energy and angular momentum, respectively. The deficit in energy and angular 
momentum is carried away by the gravitational waves, and can in principle be measured by a distant observer — it 
is therefore not surprising to find that, in the above orbital setup, the t and Lp components of the GSF are gauge 
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invariant (within a class of physically reasonable gauges). This makes it relatively easy to test the values of these 
components in numerical calculations. Indeed, both BS and SD were able to demonstrate that their numerical results 
for the t and ip components of the local GSF were in agreement with the values inferred from the fluxes of energy 
and angular momentum in the gravitational waves (the latter values were derived by BS and SD from their respective 
numerically-calculated metric perturbations, and confirmed against values tabulated in the literature). Hence, BS 
and SD's results for the dissipative part of the GSF are already well tested, and we shall not consider them further 
in this work. 

The situation with the conservative part of the GSF is more involved. Conservative GSF effects are manifest 
in a small shift in the instantaneous values of the orbital parameters away from their unperturbed values; such 
effects alter the time dependence of the orbital phases, and may have important influence on the phasing of the 
emitted gravitational waves [3- In the case of a circular orbit in Schwarzschild which concerns us here, the 
conservative effects are entirely due to the radial (r) component of the GSF. This component is not gauge invariant, 
and its value depends on the gauge in which the associated metric perturbation is expressed. A meaningful, gauge- 
invariant description of the conservative effects requires knowledge of both the GSF and the metric perturbation 
associated with it. It is not possible to check the self-consistency of the numerical calculation of the conservative 
GSF's component using energy-momentum balance considerations as done with the dissipative components. It is 
possible to test this calculation against results from the post-Newtonian (PN) literature (as done in SD by examining 
gauge-invariant quantities; see below), but such tests are only applicable for orbits with a sufficiently large radius. The 
two independent analyses by BS and SD now provide us with an opportunity to test the calculation of the conservative 
GSF effects in the strong-field regime. 

For our comparison we shall consider the two physically observable gauge- invariant quantities identified by Detweiler 
in 0, |Tl|. (The meaning of "gauge invariance" here will be made precise below; it refers to a class of gauge trans- 
formations compatible with the helical symmetry of the perturbed spacetime.) The first quantity is the azimuthal 
orbital frequency f2 (defined with respect to Schwarzschild time t), from which one derives a "gauge-invariant" orbital 
radius R [see Eq. (|4.ip below]. The second gauge invariant is U = u*, the contravariant t component of the particle's 
four-velocity. The functional form of the relation U{R) is independent of the gauge, and thus provides a convenient 
handle with which to compare between two calculations held in different gauges. 

The relation U{R) (including finite-mass terms) has been calculated in SD and used there as a basis for comparison 
with results from PN calculations. Here we will calculate U{R) based on the Lorenz-gauge results of BS. We shall see 
that a meaningful comparison between [J7(i?)]sD and [C/(i?)]BS requires two important adjustments: Firstly, one needs 
to account for the fact that the orbits in the two analyses are formally defined in two different geometries: It is an 
accelerated trajectory on a background spacetime in BS, whereas in SD it is a geodesic in a perturbed spacetime. We 
shall see that a formal connection between the two descriptions can be established by suitably relating the proper-time 
parameters along the two orbits. 

The second necessary adjustment is more subtle: The (perturbed) quantity U was shown in SD to be gauge 
invariant, for circular orbits in Schwarzschild, within a class of "physically reasonable" gauges (we shall define this 
class more accurately in the next Section). We will show, however, that the Regge- Wheeler-gauge perturbation of SD 
and the Lorenz-gauge perturbation of BS are related by a gauge transformation which does not leave U invariant. 
We will then work out explicitly a gauge transformation which brings the BS perturbation into the same class of 
"physically reasonable" gauges as the SD perturbation, and use this to gauge-adjust the value of U (the necessary 
gauge adjustment will have a monopole component only). 

The structure of this work is as follows. In Sec.|TT]we review BS's and SD's formulations of the orbital motion and 
introduce the necessary notation for our analysis. In Sec. IIIII we establish a formal relation between [/7(i?)]sD and 
[?7(i?)]BS by performing the two adjustments mentioned above. In Sec. IIVI we describe the numerical derivation of 
U{R) within the BS analysis, and compare the numerical values (adjusted for proper time and for gauge) with those 
derived in SD. This comparison is displayed in Table I. 

Throughout this work we use standard geometrized units (with c = G = 1) and metric signature ( — \—\ — h). t, r, 6, (p 
represent the standard Schwarzschild coordinates. 

II. REVIEW AND NOTATION 

BS and SD use rather different notation. For our purpose it will be useful to introduce a unified notation, which is 
one of the aims of the following short review. Our convention will be that quantities arising from the BS analysis are 
denoted with a 'tilded' symbol (e.g., X), while their SD counterparts are left 'un-tilded': 



X : quantities arising from BS, 
X : quantities arising from SD. 



(2.1) 
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We will generally denote "background" quantities (i.e., ones obtained when the metric perturbation and all GSF 
effects are ignored) with a script '0' (as in Xq). To ensure that BS and SD correspond to identical physical setups, 
we will require that all background quantities attain the same values in both analyses: 

Xq — Xq : 'background' quantities, GSF ignored. (2.2) 



A. Orbital setup and equation of motion 



In both analyses we consider a pointlike particle of mass /i(= /t), in a circular orbit around a Schwarzschild black 

0\— dap) 



hole with mass M(= M) :$> ^. The background Schwarzschild metric is denoted gal{= (jal)- At the limit /i ^ (i.e.. 



ignoring GSF effects) the particle moves on a geodesic of gf^, with a fixed Schwarzschild radius r ~ ro{= fp). (Note 

that identifying fi = fl, M — M and rp = tq is sufficient to guarantee that the underlying physical setups in both 
BS and SD are identical.) In both analyses we set up our Schwarzschild coordinate system such that the background 
geodesic is confined to the equatorial plane, 9 = n/2. 

Now consider leading-order effects arising from the finite mass of the particle. Denote the retarded linear metric 
perturbations due to the particle by hap and hap (each oc /i), corresponding to SD and BS. The two perturbations 
are mathematically distinct, since hap is taken in the Regge- Wheeler gauge 0], whereas hap is taken in the Lorenz 
gauge [H. In BS, the (regularized) back-reaction from hap is seen as giving rise to a GSF F°', which accelerates the 
particle on the background spacetime g^*^. SD, on the other hand, does not calculate the GSF directly. Instead, he 
constructs the smooth field h^^ ('R' field), which is a particular vacuum solution of the perturbed Einstein equations 
(again, in the Regge- Wheeler gauge), with the property that the perturbed orbit is a geodesic of the total (smooth) 
perturbed metric gap = gl^^ + h^p. SD then attributes back-reaction effects to the shift in the values of the orbital 
parameters with respect to their unperturbed values. 

Let us parameterize the perturbed orbits in BS and SD by their respective proper times f and r. We need 
here two different symbols, because the two orbits are formally defined in two different spacetimes: g^^^ and gap, 
correspondingly. Let events along the (accelerated) BS orbit have Schwarzschild coordinate values 5;p(f), and events 
along the (geodesic) SD orbit have "perturbed" Schwarzschild coordinate values ajp (r) [such that, for a given physical 
point along the orbit, —x^ ~ 0(/i)]. Then define the corresponding four velocities = dXp/df and u" — dx^/dr. 

(Note and u" are formally vectors in two different geometries.) In our setup we have = vP = 0, and we impose 
the circularity conditions u*", ■S'' = and du^ / dr, dvT / df = 0. 16] A given physical event along the orbit will generally 
have x" ^ a;" and ^ m", with equalities restored only at the limit /i ^ 0. An explicit relation between the 
corresponding tilded and non-tilded perturbed orbital quantities will be established in the next section. 
With the above notation, the equations of motion in BS and SD take the respective forms 

^-fr;fu''^i^=^-iF" (2.3) 



and 

du" 
~d^ 



T%yv? = 0. (2.4) 

Here the connection coefficients r^^^"* correspond to the background metric g^^p-, and F^^ correspond to the smooth 
perturbed metric gap- Note that in BS the GSF effect is entirely accounted for by the acceleration term appearing 
on the right-hand side (RHS) of the equation of motion, whereas in SD finite-/i effects are encoded in the perturbed 
values of and F^^. 

B. Gauge-invariant conservative GSF effects 

A gauge transformation in g^^^ + of the form 

x'' ^x'" =x"-FC", (2.5) 
where is of 0{iJ,), will change the SD vacuum perturbation h^^ by an amount 

Sih% = -{U.P + ^P;a). (2.6) 
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Under a similar gauge transformation in g'f^p + hal3^ the GSF in Eq. (j2.3p will vary by an amount [lo| 



(2.7) 



where an overdot denotes covariant differentiation with respect to f, Uq is the four- velocity along the background 
geodesic, and is the background Riemann tensor, evaluated at the particle. [iTl] The orbits x'^{t) and x'^{t) 

as 5^x" — (5^ip = and, taking r and f derivatives, respectively, we 



themselves will obviously change under 
find for the corresponding four-velocities 



"d7 



(2.8) 



[Note that we may use d^"/c?T and d^"/c?f interchangeably, as the two quantities differ only by an amount of 0(/i^).] 
What quantities can we construct, in our circular-equatorial orbit case, which are independent of the gauge? 
Apart from the obvious /i and r, two familiar gauge-invariant quantities are —dut/dr and du^/dr (and their 'tilded' 
counterparts) — the rate of change of the geodesic energy and angular momentum parameters, respectively. These 
quantities vanish at /z — > 0, and at the leading order [0(/i)] are entirely due to the dissipative effect of the GSF. 

For our BS/SD comparison we shall be interested in orbital gauge-invariants which are affected by the conservative 
component of the GSF. Let us restrict attention to a class of gauge transformations whose members, denoted have 
= as well as the helical symmetry 



(2.9) 



where fig = Uq /uq = \/ Mjr^ is the orbital frequency of the background geodesic. This choice is motivated as follows: 
If the original metric perturbation respects the helical symmetry of the physical (black hole -f particle) configuration, 
and is also symmetric under reflection through the equatorial plane — then a gauge transformation within the family ^ 
would retain these symmetries. The condition (|2.9p implies that, along the orbit, c?^"/dr = dS," /dr = through 0{fi). 
Equation (|2.8p therefore tells us that all components (and u") are invariant under ^. Of these four components, 
only M* and exhibit non-trivial conservative 0(/i) effects, and so could serve usefully as gauge-invariants for our 



comparison. In fact, we shall follow here SD, and utilize the alternative pair U 



and 



(2.10) 



along with their 'tilded' counterparts JJ = and Cl. The physical interpretation of the gauge invariant U is somewhat 
less obvious than that of the orbital frequency fi; it is discussed in SD. Here we shall refer to U (or U) as the time 
function along the orbit. 

The orbital frequency and time function are given in BS, through 0{p,), by 



2fiM 



(2.11) 



U = Un I 



2fi 



(2.12) 



where Uq = {1 — SAI/ro) is the time function along the background geodesic. The corresponding SD quantities 
are given by 



ro(ro-3Af) g R 
4M 



(2.13) 



where /i^^ and /i^^ ^ of course, evaluated on the circular orbit. That the quantities in Eqs. (|2.1ip - (|2.14p are gauge 
invariant under f can be readily verified with an explicit calculation, using S^tq = f'' and S^Fr = — 3/iM^'"/ [tq (ro — 3A/)] 
(see Ref. [l|), along with Eq. ([2^]) . 
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III. RELATION BETWEEN BS AND SD GAUGE INVARIANTS 

We now seek to obtain an explicit relation between fl and Cl and between U and U. We achieve this in two steps. 
In the first step we map the BS trajectory onto a geodesic y(r) defined a la SD in a perturbed spacetime (but in the 
Lorenz gauge). We define the quantities (ly and Uy associated with this geodesic, and relate them to the original fl and 
U. We then discuss the (somewhat unexpected) fact that the BS and SD metric perturbations are related by a gauge 
transformation which is not within the family ^. In the second step we therefore work out a gauge transformation 
which brings the two perturbations within a transformation, and use this to establish an explicit relation between 
fly, Uy and the SD quantities fl, U. 



A. Mapping between BS and SD trajectories 

In BS, the orbit ip(f) is an (accelerated) trajectory in the Schwarzschild background gf^. However, one can 
reinterpret this orbit, in the spirit of Detweiler and Whiting as a geodesic of a smooth perturbed geometry 
gaf3 = g^^^ + h^i^, where h^^ is the R part of the BS (Lorenz-gauge) metric perturbation haf3. (The field h^f^ has 
not been constructed explicitly in Rcf. [Ij as it is not needed for calculating the GSF in the BS approach. This field 
can, in principle, be constructed following the prescription of SD but working in the Lorenz gauge.) This geodesic is 
physically identical with the SD geodesic (r) , since h^p and h^^ represent the same physical perturbation (merely 
expressed in different gauges). We may therefore parameterize it by the SD proper time r. Let us then denote this 
geodesic by Vpir), with associated four- velocity ^"(r) = dyp/dr. Here the 'tildes' are meant to remind us that these 
quantities are described in the Lorenz gauge of BS, but one should bear in mind that both j/p and v" are defined in 
the perturbed geometry gap. 

Our goal now will be to relate the vector v°'{t) to the BS four- velocity {i"(f). This task is somewhat delicate, 
because the two vectors are defined in two different spacetimes, and so one first needs to establish a concrete mapping 
between the two trajectories Vpir) and ip(f). The mapping procedure (and its intimate relation with the gauge 
freedom of the GSF) has been discussed in detail by Barack and Ori [l^l ■ The following derivation is largely inspired 
by that work. 

We begin by obtaining a relation between the proper-time parameters t and f for a given physical point along the 
orbit. The r-interval along the geodesic ^p (t) in gap satisfies 

dr' = -{gi:^/,+h^p)dy^dfp. (3.1) 
The four- velocity v"{t) along this geodesic satisfies 

dv°' 

— +f^^(yp)z}^«^ = 0, (3.2) 

where the connections F^^ are those associated with the full metric gap, and are here evaluated at Vpir). We now 
wish to think of gf^ and h^f^ as separate tensor fields in the geometry gap. In general, of course, the splitting of 

gap into a "background" field gf^ and a "perturbation" field h^^ depends on the choice of gauge (indeed, it is the 
origin of the gauge freedom). Here, however, we fix both the background coordinates (Schwarzschild) and the gauge 
(Lorenz), and so both g^^^ and h^^ are defined unambiguously. Let T'^^'^ be the connections associated with the field 
g'^°^, and let AT'^^ = f^^ - f^^°'. Then Eq. (H^]) can be written in the form 

^ + rf'\y,)v''i'' = -Af^^(yp)^'3c^. (3.3) 
Next, define a new parameter f along the geodesic j/p (t) using 

dr'^-g^pdy^dy^ (3.4) 
[with the requirement that the "zero point" of f is chosen such that f — r = 0(/^)]; and define the new tangent vector 

dy^ dr 

u'^ ^ JpL ^i<^^_ (3.5) 
ar dr 
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We use here the symbols f and m" (in a shght abuse of notation) as we are soon to interpret these as the BS proper 
time and four velocity. For now, however, recall that f and u" are defined along the geodesic in the perturbed 

spacetime (ja/s (not in gf^); but note f and u°' are not proper time and four- velocity along this geodesic. Substituting 
v" = {df/dT)u°' and d/dr — {df/dT)d/df in Eq. (|3.3p . we now obtain 



^ + f;f{yp)uf'u^ = -Af^^(yp)w'3£t7 _ (3.6) 



where (3 = {dr/dr) (d f /dr ). The expression on the left-hand side here can be interpreted as the acceleration of a 
trajectory yp{f) (with proper time f and four-velocity -u") in a spacetime ff^^*!]- This acceleration is orthogonal to u" 
in the spacetime gf^ [by virtue of g^^^u°'u^ — —1, which follows from Eq. (|3.4p ]. so formally projecting both sides of 
Eq. (|3.6p orthogonally to u" (in the metric gf^ ) finally gives 

^ + f;f (yp)«^^^ - -{5",+u-u-^g^°M^';,,{yp)u^u''. (3.7) 

Now, Detweiler and Whiting have shown that the expression on the RHS of Eq. (|3.7p (where, recall, ^^p^ is 
derived from the R-field h^^^) is precisely the self- acceleration fj,~^F°', if one interprets this quantity as a vector in 

5a/3- Comparing then the forms of Eqs. (|3.7p and (|2.3p . we arrive at the following conclusion: If each point along 
the geodesic yp{T) is associated with a point along the BS trajectory having the same coordinate value (i.e. ip — jjp), 
then the quantities f and u" defined in the perturbed spacetime g^f^ in Eqs. (|3.4p and (j3.5p can be interpreted — and 
would have the same values — as the corresponding BS quantities f and m" = dx^/df. 

The main practical outcome from the above discussion are formulas relating the BS quantities f and to their 
counterparts r and v" defined a la SD in the perturbed spacetime. From Eqs. p.ip and p.4p we obtain 

^ = l + H (3.8) 
dr 

[through 0{p)], where 

H ^ h^0ixp)u'^u^. (3.9) 

Equation p.Sp describes the relation between the BS and SD proper-time parameters, assuming the "same-coordinate- 
value" mapping between the respective trajectories. The relation between the four-velocities then follows immediately 
from Eq. J^: 



+ (3.10) 



[again, through 0{^)]. 

Let us finally define the frequency Qy = v'^ /v* and time function Uy = associated with the trajectory ?;"(t), and 
relate these to the BS quantities Ct and U. For the frequency we have, using Eq. p.lOp . 

= = U^/U* = (3.11) 

For the time function, Eq. p.lOp immediately gives 

Uy= (i + h)u. (3.12) 



B. Adjusting the gauge 

If the BS perturbation h^^ and the SD perturbation h^^ were related through a gauge transformation within the 

family ^, then the quantities fly and Uy would have to be equal to their SD counterparts and U , since, recall, fl 
and U are gauge invariant under ^. As we now discuss, it turns out that this is not the case: The two perturbations 
differ by a gauge transformation not belonging to the family ^. 
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To see this, it suffices to examine the asymptotic form of the tt components of the BS and SD metric perturbations 
at r — > cx). In SD a gauge is chosen (within the class of Regge- Wheeler gauges) such that the perturbed metric is 
asymptotically Minkowskian: 

(r^^). (3.13) 
This is impossible to achieve with the BS Lorenz-gauge perturbation, which has a monople contribution that fails 



to vanish at r —^ oo. This monopole contribution was first derived by Detweiler and Poisson |14| . who showed that 
it uniquely describes the correct mass perturbation due to the particle: Any other Lorenz-gauge monopole solution 
would either diverge at the event horizon or at infinity, and/or fail to describe the correct mass perturbation. The 
field inherits the large-r asymptotic form of the full perturbation ha/s, which is dominated by the above monopole 
term. The perturbed BS metric has the asymptotic form 

g[',^ +hf,-^-l-2a (r^^), (3.14) 

where 

a= , ^ =, (3.15) 
Vro(ro-3Af) 

and where the term —la comes entirely from the monopole perturbation. (The other tensorial components of the 
BS metric attain their Minkowski values at r — > oo, just like the SD metric.) The difference between the asymptotic 
forms in Eqs. (I3.13|) and (I3.14|) must be accounted for by a monopole gauge transformation. From Eq. (12. 6|) we find 
that the generator of this gauge transformation must satisfy, asymptotically, ^t^t = a (as well as t,t,if = 0, since a 
monople transformation cannot depend on 95), which necessarily violates the condition of Eq. (|2.9p . Hence, the BS 
and SD perturbations differ by a gauge transformation outside the family ^. 

The above conclusion by no means suggests that either of BS/SD's metric perturbations violates the helical sym- 
metry of the physical spacetimc (they both respect it, in fact). Rather, it calls for a more careful inspection of the 
class of gauge transformations which maintain helical symmetry. Consider a gauge transformation = citS", where 
a{(x fi) is a constant. This transformation affects only the tt component of the metric perturbation, which shifts by an 
amount S^h^f. = 2(1 — 2M/r)a. If the original perturbation is helically symmetric [i.e., satisfying {dt + ^od^)h^p = 0], 
then so will be the gauge-transformed perturbation, even though the generator ^" itself does not respect the helical 
symmetry of the geometry. This suggests a more general class of "physically reasonable" gauge transformations, with 
generators given by 

C=C+atS^, (3.16) 

where ^" is any vector satisfying Eq. (|2.9p . and a is any constant (oc /i). While the generators ^ themselves are 
generally not helically symmetric, they do not interfere with the helical symmetry of the metric perturbation they act 
upon. It is not difficult to convince oneself that the class ^ is the most general class with this property. 18] 

The gauge transformation between the BS and SD perturbations belongs to the class ^, with a — a ^ 0. Crucially 
for our analysis, the two quantities f2 and U, which arc invariant under ^, are not invariant under ^ (for any a ^ 0). 
More precisely, we have (for a — a) 

S^n = -ano, S^U = aUo, (3.17) 

through 0{^). These expressions are easily derived using Eq. (|2.8p , noting d^^/dr — and hence d^^/dr = a US" 
along the worldline. Applying the gauge transformation ^" = at 5", the quantities fty and Uy of Eqs. p. lip and 
p.l2p transform as 

ily — > D,y — ai}f) = r2 — ailo 

= n{l-a)+0{ti'^), (3.18) 



and 



Uy Uy+aUo = {l + H)U + aUo 

^U{l + a + H) +0{p^). (3.19) 
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Since this gauge transformation brings the BS and SD perturbations to within a ^ transformation from one other, the 
RHS expressions in Eqs. (|3.18p and (|3.19p must be equal to the SD quantities H. and U, respectively (as, recall, the 
frequency and time function are invariant under ^ transformations). We hence arrive at the desired relations 

17 = 0(1 -a), (3.20) 



U^u(l + a + Hy (3.21) 

holding through 0{ii). 

Equations ()3.20|) and (j3.21|) explicitly relate between the BS and SD values of the perturbed frequency and time 
function. These relations do not involve the GSF, but they do require knowledge of the quantity H = h^^u^u^ /2, 

which is to be constructed from the BS perturbation. A physical interpretation of the quantity H is suggested from 
Eq. p.Sp : It describes how proper-time intervals relate to each other in the different orbital representations of BS and 
SD. It is straightforward to check that H is invariant, for circular orbits, within the class of gauge transformations ^ 
(though not within ^ for a ^ 0). 



IV. COMPARISON OF NUMERICAL RESULTS 



A. Gauge-invariant comparison formula 



The relations (|3.20|) and (|3.2ip . as they stand, do not quite yet offer a practical means by which to test the BS/SD 
numerical results against each other. The reason is as follows: The equalities expressed in these relations hold, through 
0(/i), for a given physical orbit with (perturbed) SD radius Tp and BS radius fp. In these equalities, Q and U are to 
be evaluated at Tp, while f2 and U are to be evaluated at fp. Alas, the relation between Vp and fp is not known to us 
at 0{ii): It depends on the precise gauge transformation between the BS and SD perturbations, which we have not 
solved for. Without knowledge of how relates to fp at 0{fi), it is not possible to extract the 0(/i) parts of Eqs. 
p.20p and (|3.2ip . Stated differently, we have the following problem: While the quantities and U are gauge-invariant 
(under ^), the finite-mass differences J7 — fio and U ~ Uq (which are the quantities whose numerical values we wish 
to test) are not gauge invariant, since Qq and Uq depend on the gauge (through tq). We should instead be looking at 
finite-/! corrections which are gauge independent. 

A standard solution to this problem is to express one gauge invariant in terms of the other, and, following SD, we 
shall pursue this direction here. SD introduced the gauge-invariant radius 

R={M/n^y^'^ (4.1) 

(denoted Rn in Ref. [2,]), and then expressed U in terms of R. The difference 

AU{R) = U{R) - (1 - 3A//i?)-i/2 (4.2) 

[which is ofO{fi), since (1 — 3M/i?)^^/^ — Uo + O{ii)] is then a genuinely gauge-invariant measure of the conservative 
finite-/i effect. SD derived the relation AU (R) numerically, and utilized it for comparison with results from PN theory. 
Here we shall reconstruct the relation A?7(i?) from BS quantities, and use it to compare with SD. 

Using Eq. (|2.12p to substitute for U in Eq. (|3.2ip . and then substituting the result for U in Eq. (|4.2p . we obtain 

At/ = (1 - 3Af/ro)-i/2 (^i + a + H- ^Fr^ - (1 - 3M/R)-^/^ + 0{fi^). (4.3) 

Our goal is to express the RHS here entirely in terms of the gauge-invariant radius R. For this, we need first to 
express tq in terms of R up through 0{fj,). Starting from the definition (|4.ip . then using Eqs. (|3.20p and (|2.1ip in 
succession, and finally substituting Qq^ = rf^/M, we obtain 



R^ = MQ-^ = Am-^{1 - a)-2 

„3 



'^0 



0{^?). (4.4) 
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Hence, 



ro = R 



2 R(R - 3M) ~ 



O(m'), (4.5) 



where the quantities a and Fr/fi, which are aheady 0{ii), are evaluated at R. Substituting for ro(-R) in Eq. (j4.3p and 
expanding through 0{ii), we obtain 

Af/ = (1 - 3Af/i?)-i/^ a + if) + O(M^). (4.6) 

Finally, we note that on the RHS we can replace i? ^ rg without affecting AU at leading order [since tq — R ^ 0{pi), 
and a and H are already 0{fi)]. This allows us to recast Eq. (14. 6p in a more practical form: 

At/^(l-3MAo)-/^(^o + #), (4.7) 

with corrections of 0{iJ?). Recall a = n[ro{ro — 3Af)]~^/^, and H is the quantity constructed from the BS metric 
perturbation through Eq. p.9p . 

It is interesting to point out that our final expression, Eq. (|4.7p . bears no direct reference to the GSF. This suggests 
that the 0(/i) coordinate change r — > i?, in fact, amounts to a gauge transformation which sets the conservative piece 
of the BS GSF to zero. This is further suggested by noticing Uo{ro{R)) — U{ro R) (up to terms oc a) — which one 
can readily verify starting with Uq — {1 — SM/ro)^^^^ and using Eqs. (|4.5p and (|2.12p . 



B. Numerical results 



We have reconstructed the relation AU{ro) numerically, using the BS Lorenz-gauge code, based on Eq. (14. 7p . The 
calculation of the full (retarded) Lorenz-gauge perturbation is describe in detail in Ref. In short, the method 
relies on the perturbation formalism of Barack and Lousto [isj . in which the Lorenz-gauge perturbation equations 
are decomposed into tensor harmonics, resulting in a set of 10 partial differential equations (in t and r) for each I, m- 
harmonic of the perturbation. These equations couple between different tensorial components of the perturbation, 
but are conveniently written in a form where no coupling occurs in the principal part. In BS these equations are then 
solved in the time domain (for each given /, m) using finite differentiation on a characteristic mesh. The monopole 
and dipole modes (/ = 0, 1) are calculated separately, based essentially on the semi-analytical results of Detweiler and 
Poisson 14]. The full perturbation is finally obtained by summing over a sufficiently large number of modes 1,111. 

To construct the quantity H in Eq. (j4.7p next requires us to obtain the R-part of our numerical solutions, evaluated 
at the particle. Conveniently, it is not strictly the R-part h^^ that we need, but rather the contracted quantity 
h^pU°'u^ . The construction of this quantity, through mode-by-mode regularization, is prescribed by SD in Sec. IV. B 
of [2l for any gauge related to the SD gauge through a ^-type transformation, and we employ it here in a direct 
manner. (The regular gauge transformation ^" = at5f does not affect the singular part of the perturbation, hence 
the regularization procedure described in ^ is applicable within the broader ^ family.) We construct H for a series 
of orbital radii rg, and then use Eq. (|4.7p to obtain AU for these radii. We tabulate the values thus obtained in Table 
I, alongside the corresponding SD values from Ref. The fractional difference between the BS and SD results is 
found to be similar in magnitude to the estimated fractional numerical error in the BS data. (How this numerical 
error is estimated is described in detail in Ref. [l|; the numerical error in the corresponding SD data is in all cases 
much smaller.) We conclude that the two calculations are in agreement with each other. 



V. DISCUSSION 



The agreement established here between the numerical results of BS and SD not only provides an important 
validation test for both analyses, but it also illustrates (and confirms) a few fundamental results from GSF theory. 
The following list highlights these results, (i) The GSF, as defined and calculated by Mino, Sasaki and Tanaka [5] and 
Quinn and Wald Q, can also be derived from the Detweiler- Whiting R-field 0] using the formula on the RHS 
of Eq. (|3.7p . (ii) The GSF can be defined and calculated through a "same-coordinate- value" mapping of the orbit 
from the physical perturbed spacetime onto a background spacetime; relaxing the "same-coordinate- value" rule gives 
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ro/M 


H 


Af/ (from 


BS) 




MJ (from SD) 


Rel. diflf. 


6.0 


-0.523602 


-0.296040244 


[7e - 


05] 


-0.2960275 


4e 




05 


6.2 


-0.493483 


-0.276743089 


[6e - 


05] 


-0.2767327 


4e 




05 


6.4 


-0.466941 


-0.260014908 


[6e - 


05] 


-0.2600063 


3e 




05 


6.6 


-0.443349 


—0.245359714 


[6e - 


05] 


—0.2453525 


3e 




05 


6.8 


-0.422222 


-0.232402084 


[6e - 


05] 


-0.2323959 


3e 




05 


7.0 


-0.403177 


-0.220852781 


[5e - 


05] 


-0.2208475 


2e 




05 


7.2 


-0.385906 


—0.210485427 


[5e - 


05] 


—0.2104809 


2e 




05 


7.4 


-0.370163 


-0.201120318 


[5e - 


05] 


-0.2011164 


2e 




05 


7.6 


-0.355745 


-0.192612971 


[5e - 


05] 


-0.1926095 


2e 




05 


7.8 


-0.342483 


—0.184845874 


[4e - 


05] 


—0.1848428 


2e 




05 


8.0 


-0.330239 


—0.177722443 


[4e - 


05] 


—0.1777197 


2e 




05 


9.0 


-0.280717 


-0.149362192 


[3e - 


05] 


-0.1493606 


le 




05 


10.0 


-0.244630 


-0.129123253 


[2e - 


06] 


-0.1291222 


8e 




06 


11.0 


-0.217039 


—0.113875315 


[le - 


06] 


—0.1138747 


5e 




06 


12.0 


-0.195196 


—0.101936046 


[le - 


06] 


—0.1019355 


5e 




06 


13.0 


-0.177441 


-0.092313661 


[le - 


06] 


-0.09231331 


4e 




06 


14.0 


-0.162705 


—0.084382221 


[le - 


06] 


—0.08438195 


3e 




06 


15.0 


-0.150267 


—0.077725527 


[le - 


06] 


—0.07772532 


3e 




06 


16.0 


-0.139621 


-0.072055223 


[le - 


06] 


-0.07205505 


2e 




06 


18.0 


-0.122337 


-0.062902026 


[le - 


06] 


-0.06290189 


2e 




06 


20.0 


-0.108893 


—0.055827795 


[6e - 


071 

VJ ( J 


—0.05582771 


2e 




06 


25.0 


-0.085479 


-0.043599881 


[3e - 


071 


-0.04359984 


9e 




07 


30.0 


-0.070380 


-0.035778334 


[5e - 


071 


-0.03577831 


7e 




07 


40.0 


-0.052029 


—0.026339690 


[3e - 


071 

VJ ( J 


—0.02633967 


7e 




07 


50.0 


-0.041277 


-0.020844661 


[2e - 


071 


-0.02084465 


5e 




07 


60.0 


-0.034211 


-0.017247596 


[le - 


06] 


-0.01724759 


3e 




07 


70.0 


-0.029211 


—0.014709648 


[7e - 


071 


—0.01470964 


5e 




07 


80.0 


-0.025487 


—0.012822962 


[6e - 


071 


—0.01282296 


2e 




07 


90.0 


-0.022605 


-0.011365317 


[5e - 


071 


-0.01136531 


6e 




07 


100.0 


-0.020309 


-0.010205285 


[4e - 


07] 


-0.01020528 


5e 




07 


120.0 


-0.016880 


-0.008475253 


[2c - 


07] 


-0.008475251 


2e 




07 


150.0 


-0.013469 


-0.006757093 


[3e - 


07] 


-0.006757092 


2e 




07 


200.0 


-0.010076 


-0.005050643 


[3e - 


07] 


-0.005050642 


2e 




07 



TABLE I: Comparison between numerical results from BS (Barack and Sago, Ref. l]) and SD (Detweiler, Ref. 01). The first 
column gives the radius of the circular orbit in terms of the Schwarzschild standard areal radial coordinate r [or, equivalently 
at the relevant order, in terms of the gauge-invariant radius R — see the discussion around Eqs. (|4.6|) and (|4.7|) ]. The second 
column displays the numerical value of the perturbation quantity H [defined in Eq. (|3.9|l ]. as calculated in this paper using the 
BS Lorenz-gauge code. In the third column we tabulate the numerical values of the gauge-invariant 0(/x) quantity A(7 [defined 
in Eq. (14. 2[) ]. as calculated here based on Eq. (|4.7|) and using the BS metric perturbation. The numbers in square brackets are 
the estimated fractional errors in the BS data. The fourth column displays the values of Af/ as derived using the SD code, to 
be compared with the BS values in the third column. In the fifth column we indicate the relative difference between the BS 
and SD results. This difference is found to be comparable in magnitude to the estimated numerical error in the BS data (which 
is larger than that in the SD data), providing a reassuring validation test for both analyses. 

rise to the gauge ambiguity in the GSF [l^j. (iii) A suitable gauge transformation can be made (here r — > i?) which 
nullifies the GSF (here, the conservative piece thereof). However, the information about the physical finite-/^ effect 
can then still be retrieved in full from the metric perturbation in the new gauge — cf. Eq. (|4.7p . This is a particular 
example of the general statement that the full information about the finite-/! effect is contained in the combination 
of both the GSF and the metric perturbation [To| . 

As an aside in this work, we discussed what one may mean by a "physically reasonable" gauge transformation in the 
context of circular orbits. Sensibly, a "physically-reasonable" gauge is one in which the metric perturbation admits 
the approximate helical symmetry of the black hole -|- particle system, i.e., it satisfies {dt + ^odci,)ha0 = 0. However, 



11 



it seems unnecessary to require that the gauge transformation generators 5" connecting any two such physically- 
reasonable gauges be themselves helically symmetric. If fact, the general class of such generators, denoted here ^, 
includes members ^(a ^ 0) which are not helically symmetric. The gauge transformation between the BS and SD 
perturbations — both of which being "physically-reasonable" in the above sense — is indeed generated by a vector f 
which is not helically symmetric. 

We anticipate that comparisons similar to the one discussed in this work will allow robust tests of self-force calcu- 
lations for other orbits and other spacetimes (e.g., Kerr) when such calculations are available. The essential elements 
of our formal discussion are directly applicable to other orbits and geometries. Most important, Eqs. p.Sp and (|3.10|1 . 
which describe the mapping of the orbital elements from a BS-type background spacetime to an SD-type perturbed 
spacetime, hold quite generally for any orbit in any black hole spacetime, and could form a basis for future compar- 
isons. The major challenge in any such future comparison would remain to devise a suitable set of gauge-invariant 
quantities. 
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